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Transport through coupled quantum dots in a phonon bath is studied using the recently developed 
real-time renormalization-group method. Thereby, the problem can be treated beyond perturbation 
theory regarding the complete interaction. A reliable solution for the stationary tunnel current is 
obtained for the case of moderately strong couplings of the dots to the leads and to the phonon 
bath. Any other parameter is arbitrary, and the complete electron-phonon interaction is taken 
into account. Experimental results are quantitatively reproduced by taking into account a finite 
extension of the wavefunctions within the dots. Its dependence on the energy difference between 
the dots is derived. 
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Introduction. Today quantum dot systems allow a de- 
tailed study of many physical phenomena, like Coulomb 
blockade , Kondo effects [[| [| , interference effects [Q . 
As in these structures quantum states can be manipu- 
lated, they also may have an application in future quan- 
tum gates ||. Quantum dot systems can typically be 
characterized by only a few parameters, which are ex- 
perimentally controllable. Theoretically, these systems 
may be described by basic models, which capture the es- 
sential physics and can be investigated using standard 
methods of many-particle theory. However, if one deals 
with out-of-equilibrium situations and strong coupling, a 
theoretical analysis becomes difficult. 

Such an out-of-equilibrium problem was recently stud- 
ied, when the stationary tunnel current through a double 
quantum dot in a phonon bath was measur ed @. There 
the influence of the phonon environment was examined 
at low temperature. The thermal energy of the environ- 
ment is always a source of unwanted transitions in quan- 
tum dot devices. Even at zero temperature spontaneous 
emission of phonons gives rise to inelastic transitions, 
i.e. they occur between dot states of nonequal energy. 
In the experiment in Ref. [Q the inelastic contribution 
to the tunnel current through the double dot was stud- 
ied. A first theoretical interpretation of the experimental 
results focused on the interaction of the dots with the 
phonons, which is analogous to the spin-boson model ||. 
They accounted for the coupling to the leads perturba- 
tively and studied the electron-phonon problem using an 
approximation, which corresponds to the noninteracting 
blip approximation (NIB A) for the spin-boson model . 
Thereby a better understanding of the shape of the in- 
elastic current spectrum was achieved. It was shown that 
the interference of phonons interacting with the electron 
densities at the two dots leads to an oscillating structure 
in the current spectrum. However, a quantitative com- 
parison with the experiment has not yet been possible. 
Especially, the unexpectedly large inelastic current of the 
experiment could not be explained. 

In this paper we present a reliable solution for the sta- 



tionary state of the double dot system. Since we deal with 
a nonequilibrium situation, where both the coupling to 
the phonon reservoir and the coupling to the leads have 
to be considered as strong, we use the recently developed 
real-time renormalization-group (RTRG) method. It 
has successfully been applied to both equilibrium |[l 12 
and nonequilibrium problems |]l3| , |l4f , including the spin- 
boson model Jl4|,[l5| ■ Applying this approach to the dou- 
ble dot system both the electron and the phonon reser- 
voirs are integrated out in a RG procedure, i.e. beyond 
perturbation theory. Thereby, the external voltage is ac- 
counted for properly, and the level broadening induced 
by the coupling to the leads is included in this method. 
Thus, we can also consider the case, where the energy 
difference between the dots is of the same order as the 
external voltage. Furthermore, we do not have to in- 
clude an additional cutoff parameter, which simulates 
the level broadening generated by the leads. Moreover, 
this formalism may treat any form of dot-phonon inter- 
action. Therefore, we are able to account for the full 
electron-phonon interaction, i.e. including interaction 
terms, which involve a tunneling between the dots ("off- 
diagonal interaction terms"). In the quantitative anal- 
ysis the offdiagonal interaction terms lead to a strong 
dependence of the current on the extension of the wave- 
functions within one dot. We find that the variation of 
this extension with the energy difference e between the 
dot levels has to be accounted for. By fitting the result 
for the current with the experimental data, we obtain the 
width of the electron density within one dot as a function 
of e. Furthermore, for the physically relevant situation, 
where the coupling of the dots to the leads and to the 
phonons is only moderately large, the RTRG approach 
does not deal with any parameter restriction. This is in 
contrast to the NIBA, which is valid only for sufficiently 
high temperature. 

Model. Our model consists of two coupled quantum 
dots (I and r, respectively). Each dot is coupled to an 
electron reservoir with the chemical potentials /xj and fi r , 
see Fig. [l| We consider the case realized in the experi- 
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ment |7) , where the external voltage V = pi — p r is much 
smaller than the Coulomb charging energy U. Thus, due 
to Coulomb blockade the double dot cannot be charged 
with more than one additional electron. In the experi- 
ment a strong magnetic field was applied perpendic- 
ular to the dots. Thus, we assume spin polarization here 
and omit the spin index. We denote the many-particle 
ground states, where an additional electron is in the left 
(right) dot, by \l) (\r)) and neglect any excited states. 
Therefore, together with the uncharged ground state |0), 
there are only three possible states of the double dot. 
The total Hamiltonian H for the system can be written 
as a sum of the dot Hamiltonian, the contributions of the 
electron reservoirs and the phonon bath, and the interac- 
tion parts stemming from the coupling to the leads and 
the electron-phonon interaction: 

H = Hd + H res + H p h + H e - res + H e - P h . (1) 

The dot Hamiltonian Hd reads 

ffd = e I |l)</| + e r |r>(r|+T c (|l)<r| + |r)(Z|) , (2) 

where e; (e r ) are the ground state energies of | Z) and 
the coupling between the dots is described by the tunnel 
matrix element T c . The reservoir contributions are given 
by 



k k 

H ph = 2J u q a\a q . 



(3) 
(4) 



Here, the operators ct (c k ) create (annihilate) an electron 
with the energy e k in the left lead, whereas the creation 
(annihilation) operators d k (d k ) refer to the right elec- 
tron reservoir. Analogously, a q (a q ) create (annihilate) 
a phonon with the wavevector q and the frequency uj q . 
Here and throughout the paper we set h — 1 . The double 
dot is coupled to the external leads by the parameters V k 
and Wk- 



H e . res = (VkC k \l){0\ + V k *\0)(l\c^ 

k 

+ Y / (Wkdk\r}(0\+Wm(r\d{) ■ (5) 



The electron-phonon interaction consists of a diagonal 
part, which is characterized by the coupling constants a q 
and /3 q , and an offdiagonal contribution with the param- 
eter 7 g : 

He- P h = 52{a 9 \l)(l\+0 q \r)(r\){al + a„ q ) 
q 

+5> 9 (IO(H + k></|)(4 + (6) 



The above interaction coefficients are given by JsJ 



a q = X q {l\e^ s \l), 
(3 q = \ g (r\e i * s \r), 
lq = X q (l\e^\r) , 



(7) 
(8) 
(9) 



where X q is the matrix element for the interaction of 
2DEG electrons and phonons. The phonons are assumed 
to be three-dimensional acoustical phonons Q|. It then 
follows for the interaction m 



\K\ 2 = 9 



V\q\ 



and the dispersion reads 



c s |<z1 



(10) 



(11) 



Here, we introduced c s as the speed of sound in the 
medium, V as the volume of the crystal, and the di- 
mensionless coupling constant g For the evaluation 
of Eqs. (0) - (^) we model the electron densities pi(x) 
(p r (x)) within one dot by Gaussians, which are peaked 
around the dot positions xi [x r = xi + d) with a width 
\Ax\ = y/S/2a: 



Pi(r)(x) = — J 



1 



.3/2 



^2 — 



(12) 



The finite width a leads to a high-energy cutoff D = 
c s /(j for the coefficients a q , (3 q and j q . We include this 
cutoff in an exponential form, so that we end up with the 
following interaction coefficients: 



uc q 



X q e lqx 'e- 



lq = X q e lq{ — ^e- 



(13) 
(14) 
(15) 



A simple form of the Hamiltonian, which shows the anal- 
ogy with the spin-boson model, is obtained by shifting 
the bosonic field operators. One introduces the unitary 
transformation 



U — exp 



E 



2uj 



h t 
— -al 



Pa 



so that 



Ua a U^ = a a 



2uj 



2u>„ 



Thus, our final Hamiltonian H = UHU^ reads 

H = Hq + Hb + Hy , 

\m\-V){r\) 



(16) 



(17) 



(18) 



Ho 
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+T c cff (|0(r| + |r)(Z|) 

+£(|Z)<Z| + |r)(r|-|0><0|), (19) 
H B = ^ 6 k C k C k + X! 6kd k dk + X! W 9 a « a 9 ' ( 2 °) 



R v = : mi* 



(21) 



where we have used Eqs. (|ll| ) and (|13|) - (|15[). Further- 
more, for simplicity we have set (e; + e r )/2 = and have 
introduced the parameters 



e; - e-r , 



T c =T C - 2gu d e arctan - — , 

2u>d 

9 ( D 
E = — — D + ujd arctan — 
4 V u d 



and 



\d\ 



(22) 
(23) 

(24) 



(25) 



Thus, the tunnel amplitude T c has to be replaced by 
a smaller effective T° ff , which is due to the off diagonal 
electron-phonon interaction. One already recognizes that 
the reduction of T c strongly depends on the width of the 
electron densities a — \d\oJd/D. Finally, in view of the 
RTRG method, we have written the interaction part Hy 
as normal ordered products of local (dot) operators g^ 
and environmental operators j M . They are defined by 
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9b 2 

jb-2 



(|f>(l|-|r)<r|) 



= (IOvi + |r)(Z|) , 



(26) 
(27) 

(28) 
(29) 



neglecting the offdiagonal electron-phonon (62) interac- 
tion. 

RTRG approach. The RTRG approach is based on the 
formally exact kinetic equations 



</)(i) = Tr 



dt' ^i{t-t')p{t') 



p(t) + iL p(t) = / dt' E(f - t')p(t') 
Jo 



(35) 
(36) 



where (I)(t) denotes the time-dependent expectation 
value of the current through the double dot, and p(t) 
is the time-dependent reduced density matrix, which is 
the trace over the bath degrees of freedom Tr# of the 
full density matrix. In contrast, Tro denotes the trace 
over the local (dot) degrees of freedom. The bath de- 
grees of freedom enter the equations via the integral 
kernels E and Ej, which are defined by sums over ir- 
reducible diagrams, see Ref. (l(| for details. Introduc- 
ing the Laplace transforms f(z) = J °° dte lzt f(t) of the 
time-dependent functions f(t) and using the identity 
limt^oo f(t) = —i lim 2 ^o z f( z ) leads to the following ex- 
pression for the stationary tunnel current I s t'. 



I st = Tr [E/(z = 0)p st ] 



(37) 



where the stationary reduced density matrix p st is deter- 
mined by 



(L o + «E(z = 0))p sf = 0. 



(38) 



In the following we outline only the main steps of the (i) 
derivation of Eqs. (§|) and @, and of the (ii) RTRG 
technique, by which E and E/ are calculated. 

(i) The quantity of interest is the current through the 
double dot. It is given by the expectation value of the 
operator 



I = ie 



J2{v:\o)(U-v k c k \i)(o\ 



(39) 



9b 3 
jb 3 



-5IOXOI 



( a 9 + a?) (4 + a -<< 



9+l=9U = \l)(0\ 



9+r = 9-r 



\r)(0\ 



Ck 



J+r = 3- 



J2 w * d k 



(30) 
(31) 

(32) 

(33) 
(34) 



Therefore, the interaction index fx runs over the bosonic 
indices 61,62,^3 and the fermionic ones +1, — I, +r, — r. 
From Eqs. (|l|) - @ the spin-boson model is recovered 
by omitting the electron reservoirs and the interaction 
with them, excluding the state |0) and accounting only 
for the bosonic interaction b\. The latter corresponds to 



with e being the elementary charge. For the calculation 
of (I) (t) we introduce the Liouvillian L — [H, ■] _ and the 
superoperator A] = •] + , where [•, ■]_(+) denotes the 
(anti)commutator. The expectation value of the current 
can then be written as 



(I){t) = -*Tr [Aje^p^pl 



(40) 



Here, we assumed a factorized density matrix at t = 0, 
where p(0) is the initial dot density matrix and p e ^j is 
the equilibrium distribution of the electron reservoirs and 
the phonon bath. To evaluate the above expression the 
propagator exp (—iLt) is expanded in the interaction part 
Ly = [Hy, •]. According to Eq. ( fjll) this can be written 
as 



Ly = £ : GIJV 



(41) 



pi' 
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where the index p = ± denotes, whether the interaction 
takes place on the forward or backward propagator, i.e. 



G- = • (-<?„) 



(42) 
(43) 



Correspondingly, we write ^/ = Ep : (^;+i^+i + ^/-!^ ? i) : 
with 



= ||0<o|- , A7 +I = .||I><0|, 



.4 



7-i 



|0)(Z| 



.4 



i-l 



\0){l\ 



(44) 
(45) 



The trace over the bath degrees of freedom can then be 
performed by application of Wick's theorem. In this way 
one obtains a series of terms, where vertices G?J and A p I±l 
of the local system are connected by pair contractions 
7™' {t) = Tr B [J£J?'pb] of the bath. Denoting the sum 



over all irreducible diagrams, which contain the current 



superoperators A p ±l , by E/ then leads to Eq. (35), see 



Fig. |2j. Similarly, Eq. fl36| ) is obtained by introducing the 
object S, which is defined as the sum over all irreducible 
diagrams involving only the vertices G p . 

(ii) The objects and are calculated by a 

renormalization group procedure. Short time scales of 
7u«' W are integrated out first by introducing a short- 



time cutoff t 



into the correlation function , (t) 



pp i{t, t c ). In each renormalization-group step, the time 



scales between i 



_ and t c + dt c are integrated out, starting 
from i c = and ending at t c — 00. As a consequence, 
one generates RG equations for £(2) (£/(z)), Lq, G^, 
and the boundary vertex operators A? (Aj )and (de- 
fined as the rightmost and leftmost vertex of the kernel 
S(z) (S/(z))), for the explicit form of the RG equations 
see the Appendix. 

Within the scheme of a perturbative RG analysis, 
the generation of multiple vertex superoperators is ne- 
glected here, as in Ref. (l^jlj]. Due to the small real- 
istic value of g (g — 0.05 for GaAs) the neglecting of 
double- and higher-order vertex objects is justified and 
our approach leads to very reliable results, see Ref. p4[ , 
where we solved the spin-boson model for couplings up to 
a = g/2 < 0.1. .0.2. Furthermore, since also the coupling 
to the leads is treated nonperturbatively, the induced 
level broadening is accounted for properly. In contrast, in 
a perturbative analysis this effect has to be accounted for 
by an additional cutoff parameter. Also note that, in our 
method there is no restriction regarding the temperature, 
whereas the NIBA is only valid for the parameter regime 
(A? + e 2 ) 1/2 < T @. Here, A r = 2T C (2T C /^) Q /( 1 - Q ) 
is the renormalized tunnel amplitude of the spin-boson 
model. Finally, using the RTRG we are also able to 
account for the offdiagonal electron-phonon interaction, 
which is important for small D/uid, see Eqs. (^3|), (5: 



Results. We solve the set of ordinary differential equa- 
tions, Eqs. (|i]) - J68|), numerically. The stationary 
tunnel current then follows from Eqs. (^) and (|3^). 
Our choice of the parameters corresponds to the experi- 
ment, where a GaAs structure was used at the temper- 
ature T = 23mK = 1.98fieV with an external voltage 
V = 14Q^ieV 0. For GaAs we have c s = 5000m/s and 
g = 0.05 JXt| . The distance between the dots is estimated 
as d — 200 • 10~ 9 to B, which leads to tod — 16.5^ey. 

The result for T c = T ; = T r = lfjteV, A = D r = 
ImeV and D = 100/ieV respectively D = 150^,eV, 
which corresponds to the parameters studied in Ref. ||, 
is shown in Fig. |l The external voltage V gives rise 
to a finite stationary tunnel current through the dou- 
ble dot. The elastic current can be seen around e with 
a width depending on the coupling to the leads T; = T r 
and the internal tunnel amplitude T c . There the phonons 
do not participate in the tunnel process. Due to the cou- 
pling to the phonons there is also an inelastic current, 
where phonons are emitted (e > 0) respectively absorbed 
(e < 0) during the tunnel process. For increasing rj( r ) 
the width of the elastic current grows, while an increased 
coupling constant g leads to a larger inelastic current. 
The effect of the finite voltage V can be seen in Fig. ||, 
where for e > V the tunnel current drops to zero. Fur- 
thermore, Fig. H shows that the offdiagonal interaction 
leads to a larger inelastic current. This effect is increased 
with decreasing D, see also Eqs. ([33]) - (prf). We will see 
below (Fig. ||) that, due to Eq. ([23D, the offdiagonal in- 
teraction also suppresses the elastic current. Eventually, 
in Fig. ^ one also recognizes the oscillations stemming 
from the interference of the phonons interacting with the 
two dots ||. 

Let us now study the current quantitatively in com- 
parison with the experiment. For this it is necessary 
to choose realistic parameter values for T Cl Tir r \, D and 
Du r ) ■ From changing the bias polarity in the experiments 
the ratio T r /Ti « 0.5..1 was found M. To determine the 
couplings T c and T r , the experimental data are compared 
with the result of Stoof and Nazarov jlg| 

T 2 r,. 

Ist = r c 2 (2 + r r ./r c i ) + r2/4 + e 2 ' (46) 

which is valid for no electron-phonon interaction [g = 0). 
A good agreement of the elastic current is found for T c — 
O.UA^eV and r ; = T r = 3.5fieV, see Fig. f§ 

However, due to the absent electron-phonon interac- 
tion the influence of the finite extension of the electron 
densities within the dots is also disregarded in the Stoof- 
Nazarov result. In contrast, our method accounts for this 
extension by the high-energy cutoff D. From Eq. ( |23| ) we 
see that for a finite D the tunnel amplitude is effectively 
reduced. Therefore, the Stoof-Nazarov result underesti- 
mates the value of T c . 

From the large inelastic current in the experiment one 
can conclude, that in fact a finite value of D was re- 
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alized. It turns out that T c w 0.375 allows sensible 
fits. First, in Fig. || our results for T c = 0.375/xey, 
r ; = r r = 3.5/xeF, D t = D r = ImeV and D = 70/xeV 
respectively I? = lOO/xeV^ are compared with the experi- 
ment. One recognizes that with decreasing D the larger 
overlap of the dots' wavefunctions leads to a stronger im- 
pact of the offdiagonal electron-phonon interaction. The 
elastic current is suppressed, whereas the inelastic cur- 
rent is increased. The deviations from the experiment 
show, that there is an e dependence of the width of the 
electron densities, which we have to account for in order 
to achieve agreement. In Fig. |^ we show a fit of the width 
a = dujd/D, which is based on the experimental results 
for I st . One recognizes that for larger absolute values of 
e the electron densities are more sharply peaked. The 
asymmetry is due to the finite external voltage V. For 
e < the state \l) lies in a deep potential well, thus this 
energetic separation of the two quantum dot levels and 
the leads causes a very small overlap of the wavefunctions 
within the dots. On the other hand for e > neither dot 
level lies in a deep potential well, however, an increasing 
energetic separation e again leads to more sharply defined 
electron densities. 

In Fig. ^| one also recognizes that the structure on the 
emission side (e > 0) of the current spectrum observed in 
the experiment does not stem from interference effects of 
the phonons. In fact, the oscillations generated by this 
interference occur on a much larger energy scale than the 
structure found in the experimental curve. In contrast 
our results show that the bump on the emission side of 
the current spectrum is due to another mechanism: on 
the one hand we see in Fig. || that for small e > the 
inelastic current grows with increasing e. On the other 
hand however, the electron densities are simultaneously 
sharpened, so that the cutoff D is increased. For larger 
e this again reduces the inelastic current. 

In summary, we have applied the RTRG method to 
the coupled quantum dot system in a phonon bath in 
nonequilibrium. By accounting for both the coupling to 
the leads and the coupling to the environmental phonons 
nonperturbatively we achieved a reliable solution for the 
stationary tunnel current. For the first time both the 
elastic and the inelastic current of the experiment could 
quantitatively be reproduced. Our analysis shows the 
importance of the finite width a of the electron densities 
within one dot, and for the experiment, the dependence 
of a on the energy difference e between the dots was 
calculated. 

Acknowledgments. We acknowledge useful discussions 
with T. Brandes and K. Schonhammer. We also thank 
T. Fujisawa for the experimental data. This work was 
supported by the "Deutsche Forschungsgemcinschaft" as 
part of "SFB 345" (M.K.) and "SFB 195" (H.S.). 



Appendix 

To obtain the explicit form of the RG equations one 
has to make a choice of the t c dependence of the bath 
contractions 7^/ Since the bosonic bath contrac- 
tions lf p bk {t) (j,k = 1,2,3) correspond to those of the 
spin-boson model, we choose the t c dependence as in 
Refs. UH): 



7% bh (t,t c ) = -(R jk (t)G(t-t c ) 
+ip'S jk {t)Q{t-t c ). 
The functions Rjk{t) and Sjk{t) are defined by 

From Eqs. (|lf) - © we obtain 

R n (t)=gRe irT coth (nT(t - i/D)) 

Ud , f sinh {nT{t - l/ujd - i/D)) 
2 U \smh{irT(t+l/uj d -i/D)) 

V«3 



Suit) = glm 



-{{t-i/Df-l/^) (t-i/DY 
R12 — R21 = S12 = S21 = , 
-R13 = R31 — <Si3 = S^i = , 



R22W 
S22W 



_ 9 -0/1*4 



Re 



= _9 -D/u> d 



lm 



TrTcoth (nT(t-i/D)) 
1 



£23 0) =Rsa(t) = -2 Re 



(t- i/D) 2 



ln / sinh (nT(t-l/2u d -i/D)) \ 
X n Vsinh (wT(t + l/2co d -i/D))) 



S 23 (t) = S 32 (t) = -glm 



e -D/2uj d 

.{{t-i/D)*- If 
R 33 (t) =gRe 7rTcoth(7rT(i -i/D)) 



LUd , / sinh {nT{t - l/u d - i/D)) 
~2 11 Vsinh (irT(t+l/uj d -i/D)) 



S33 (t) = fflm 



l/^ d -2{t-i/Df 



{{t-i/DY -1/ul) {t-i/DY 



(47) 



(48) 



(49) 

(50) 

(51) 
(52) 

(53) 
(54) 

(55) 
(56) 

(57) 
(58) 



Here, we assumed that the high-energy cutoff D is much 
larger than the temperature T. In Eqs. ([33]) - ((5^) one 
again recognizes that the influence of the offdiagonal in- 
teraction strongly depends on the width of the electron 
densities a — \d\uj d /D. The fermionic contractions can 
be written as 

T?Mj,{t) = 6 n> _ n ,6 fJ , Q ( 7n/ (i) +7-„/H)) 

+ j(lr,f{t)-l-r,f{-t))^ , (59) 
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where the indices r\ and / run over ± and I, r. We intro- 
duce 



T l (e) = 2nJ2\Vk\ 2 5(e-e k ), 

k 

T r (e) = 2nY / \Wk\ 2 S{e-e k ) 1 



so that we obtain 



7*j/ (*) 



2sinh (nT(t - i/D f )) 



(60) 
(61) 

(62) 



Here, we introduced a bandwidth Df of the reservoir / 
and assumed, that Tf(e) « const, holds. The cutoff- 
dependence of the fermionic contractions is chosen as 



(63) 



The final RG equations then read (for a more detailed 
derivation we refer to Rcf. jig]) 



dtr 



dL 

dt c 



dG£ 

dt c 



= E -^(^fc(*c)Aft 6j (i o )(£ -«) 
piiP2 ,i,fe 

+P2^(i c )^) )b .(4))^ 

- E 7-%/(*c)^ lP M-^ / (t c ) J B-, (64) 

Pl,P2,V,f 

= E (^ fe (ic)[G^(t c ),L ] 

Pi,P2,j,k 

+V2S 3k (t c )Gl)(t c )) Gil 
4 E 7^% / (*c)^ lP2 G^ / (t c )G^ / , (65) 

Pl,P2,V,f 

= E ^fc(*c)(G^(*c)^-GJG^(t e ))G^ 

- 4 (G?;(t)^-^^(t) 



II 



X 



(R jk (tc)[L ,G p b l(t-t c )} 
+p 2 S jk (t c )GH(t-t c )) 

pi,P2,nJ 

v7^ PlP2 G% f {t)Gl) G%(t - t c ) , (66) 



dt c 



= E ^-*(*c)U^ 6i (t c )G? 



Pl,P2,j,k 



-A^Glld^Gll 



rff 



(^ fc (t c )[L ,GS£(t-t c )] 



+P2^- fe (f c )G^(f-f c )) 

+ E 

Pl,P2,T]J 

x jf '* K/)/ pi?2G -w(*) 

-^^^^(t)^) G*£(t - * c ) , (67) 
^= E RMGl){t c )GlBll-if°dtGll{t) 

C Pl,P2,j,k 

xGl(R 3k {t c ){L -z)BH{t-t c ) 

+P2S jk (tc)B p b 2 k (t-t c )) 
- E A%/(*c) 

Pl,P2,Vj 

x dt v ^a^G p ^ f (t)G^j(t-t c ). (68) 
Jo 

The interaction picture is defined by G^(f) = 
e ii D t G » e -iLot ) A?^ (t) = e lz *A^ )Ai e- 4io *, and = 



e lLot BPf 



The function 77^ and the superoperator 



<r pp account for additional signs arising from the com- 
mutation of fermionic field operators. They are given by 



v pp 



-pp' for \x fermionic 
1 else 



(69) 



pp' for iV s - JVy = odd 
1 for TV, - N.i = even ' 1 > 



We also note, that in the limit of large Di and D r the 
first part of the fermionic contractions in Eq. ( |59| ) can be 
written as 



^(7n/(*)+7-,/H)) = ^(*)- 



(71) 



Thus, for Di, D r — > 00, this contribution to the differen- 
tial equations Eqs. (|64|) - fl6Sj ) can be incorporated in the 
initial conditions. 
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<I>(t) = 




t' 



t 



FIG. 2. Diagrammatic expression for (I)(t). The two lines 
of the Keldysh contour are put together to one line. The 
irreducible diagrams in Ej include the leftmost vertex super- 
operator at the time point t' and the current superoperator 
Ai at the time point t. £/ acts on p(t'), which is represented 
by the horizontal line. 
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FIG. 3. Stationary tunnel current as a function of e for 
T c = Ti = T r = lfieV and D t = D r = lmeV . and 
D = 100/ieV respectively D = 150/ieV . Solid line: 
D = 150^eV. Dashed line: D = 150^eV, -y q = 0. Dotted 
line: D = 100^eV. Dot-dashed line: D = lOOpeV, j q = 0. 



a. 



FIG. 1. The double quantum dot may be charged with one 
additional electron in the left or right dot. The corresponding 
states |i) and \r) are coupled by the tunnel amplitude T c . 
The couplings to the leads are given by F;( r ). The energy 
difference between the quantum dots is e = e; — e r , and there 
is an external voltage V = — fi r . The interaction with the 
acoustical phonons (dashed lines) , consists of a diagonal part 
with the coefficients a q and /3 q , and an offdiagonal part with 
the constants j q . 



"s/PA 



Experiment 
Stoof-Nazarov 




FIG. 4. Stationary tunnel current as a function of e. Solid 
line: Experiment. Dashed line: Stoof-Nazarov result for the 
case of no electron-phonon interaction, with the parameters 
T c = 0.124/ueV and T l =F r = 3.5^eV. 
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FIG. 5. Stationary tunnel current as a function of e. Solid 
line: Experiment. Dashed line: RTRG with T c = 0.375/jeV, 
r ; = T r = 3.5fj,eV, Di = D r = lmeV and D = 70fieV. 
Dotted line: RTRG with T c = 0. 37 '5 fieV , F t = T r = 3.5/ueV, 
Di=D r = lmeV and D = KXVieV. 
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FIG. 6. The width of the electron density within one dot, 
a, as a function of e. T c = 0.375/ie^, T; = T r = 3.5/ieV and 
Di — D r — lmeV. 
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